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Extensive studies have demonstrated that finite-range rcgularization (FRR) offers sig- 
nificantly improved chiral extrapolations for lattice QCD. These studies have typically 
relied on selecting the finite-regularization scale based upon phenomcnological input. 
Here we report on a preliminary investigation of a procedure to determine a preferred 
range of FRR, scale based on nonperturbative lattice results — without any phenomeno- 
logical prejudice. 



1. Background 

Modern lattice QCD results are beyond the power-counting regime (pcr|H21 of 
the chiral expansion. Any results from implementation of the chiral effective field 
theory (EFT) in this domain are therefore model dependent. From the perspective 
of lattice practitioners, the goal is focussed on ab initio studies of nonperturbative 
QCD — without phenomenological input. Thereby, even when working outside the 
PCR, it can be attractive to choose dimensional regularisation (or a similar minimal 
subtraction-sfyZe scheme) to avoid phenomenological bias. As results are dependent 
upon this choice of how the EFT is applied, then this application is implicitly a 
model. Nevertheless, the ab initio demand is maintained by not introducing addi- 
tional information from prior knowledge. 

Finite-range regularisation (FRR)^ has been demonstrated to offer stable and 
robust chiral extrapolation of lattice simulation results performed at moderate 
quark masses. The FRR procedure necessarily involves introducing a preference 
for a regularisation scale that is "low" in favour of one that is infinite — whereby 
this assumption of prior knowledge means one can no longer assert ab initio. 

Here we highlight an initial investigative study in search of a technique which 
will enable the results of lattice simulations to select a preferred regularisation scale. 
By removing the need to artificially select the regularisation scale by hand, one can 
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utilise the benefits of FRR and maintain ab initio status. 
2. Convergence and the power-counting regime 

The the success of effective field theories, such as chiral pertubation theory, lies in 
the presence of a small expansion parameter and the reliance of a series expansion 
of natural-sized coefficients. A simple example function is the geometric series 

— — = 1 - x + x 2 - x 3 + . . . (1) 

1 + x w 

The radius of convergence of this series expansion about x = is |a;| < 1. Given 
a truncated expansion, the deviation from the exact result, at a given value of x, 
can be estimated by considering the size of the first neglected term in the series — 
without necessarily having knowledge of the full result. 

This toy example can provide insight into the natural convergence radius of the 
quark-mass expansion of the nucleon mass. This expansion can be written as 

M N = M% + c 2 ml + Xrrml + . . . , (2) 

where oc ^Jm q . The term is nonanalytic in the quark mass (ie. ~ m^ 2 ), 
and it's coefficient is know model-independently, 

3 9a ^ 

Xrr- 32?r/ 2. (6) 

Strictly, qa and f w are to be evaluated in the chiral limit. Assuming the phys- 
ical values do not differ too significantly, this suggests a numerical value of 
Xir — — 5.6GeV~ 2 . The term in m% is also known phenomenologically through 
the pion-nucleon sigma term, whereby^ 

a N = ml^- ~ 45MeV . (4) 

Using this and the physical nucleon mass, together with the leading nonanalytic 
term gives, c 2 ~ 3.5 GeV -1 and M% ~ 0.89 GeV. 

Using these numerical estimates in Eq. j2j one finds that the second and third 
term reach 100% of the leading term at pion masses of 0.51 GeV and 0.54 GeV, 
respectively. This suggests that at a pion mass of the order ~ 0.5 GeV one is 
at the radius of convergence of this series — where an infinite number of terms are 
required to reproduce the exact result. Of practical importance is how precisely can 
the curve be reproduced with a finite number of terms. 

The first few terms of the expansion Eq. |2]) looks very much like the geometric 
series, setting rriR = 0.54 GeV, 



July 2, 2009 15:39 WSPC/INSTRUCTION FILE nz'rho 



Towards selecting a finite-range regularization scale 3 

The suggests that using Eq. |1| with x = m n /mp/ (ignoring the term linear in x) 
is a good way to estimate the precision of a given truncation. Demanding that the 
series expansion in Eq. ([1]) up to x 3 is accurate to the 1% level, this limits one to 
x < 0.34. Thereby 1% precision is enabled up to m^/m^ < 0.34 or < 0.18 GeV. 
Continuing the natural-size argument means that at order m\ demanding this level 
of precision restricts one to < 0.23 GeV. Expansions to higher order in the 
nucleon mass have also recently been revisited^lSI, but even going to m% or m\ 
still only provides a valid expansion up to < 0.27 and 0.30 GeV, respectively. 
Further, going to increasing order also comes at the cost of fitting more parameters 
and thereby requiring even more simulation results below the quoted thresholds. 

3. Numerical approach to FRR scale determination 

With the aim of performing precision ab initio studies, it is evident that lattice 
results from beyond the (1%) PCR must be used. The first option would be to 
just to use the truncated forms of a DR expansion. In using lattice results up to 
m n ~ 0.5 GeV to constrain an extrapolation based upon a fourth-order expansion 
would lead to a systematic uncertainty of 35% (based on the geometric series dis- 
cussed above). A second option is to introduce a finite-range regularisation scale 
with the aim of both incorporating more data, to improve statistical precision, and 
minimising the systematic uncertainty. 

The physical interpretation for the success of FRR lies in the suppression of 
rapidly- varying chiral logs at moderate quark masses. Mathematically, the success 
is based on using a separation of scales to introduce an effective resummation of 
higher-order terms in the chiral series, while maintaining the model-independence 
of the expansion to the order one is working. 

To continue with the nucleon mass example, the chiral expansion, Eq. ([2]), can 
be rewritten in an unrenormalised form 

Mjv = a + a 2 m^ + a±m% + S(m w , A) , (6) 

where £ denotes the meson-loop corrections, which contains the chiral nonalytici- 
ties, and a, are regularisation-scale dependent expansion coefficients. In adjusting 
the regularisation scale A, one can shift strength bewteen the loop contributions 
and the residual series composed of the a^s. FRR works well because with an ap- 
propriately chosen A, the residual series typically shows much better convergence 
than the renormalised form. We wish to exploit this feature in using lattice results 
to choose an "optimal" range for A. 

Further, any criteria to determine A should produce certain desired limits. 
Firstly, if one really is using data in the PCR, then there should be no perefer- 
ence for a finite A, and A — > oo should display equally good covergence properties. 
In fact, the data should always reject a A which is too low. In the PCR, then let- 
ting A — > would give a curve inconsistent with the nonanalyticities of the data. 
Further, if the expansion is extended to quark masses beyond the range accessible 
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within FRR, then a solution consistent with A — > is a clear indication that the 
data is inconsistent with the chiral behaviour. In this domain, it may also be evident 
that no amount of tuning A can stablize the convergence. 

The technique to constrain A we investigate here, is to numerically adjust A 
such as to ensure the residual series is satisfactorily converging over the range upon 
which the lattice results are fit. In the case above, Eq. ([6]), a preferred value of A is 
one which gives numerical evidence of convergence of the series 

{a , a 2 ml maK , a 4 m^ max ,...}, (7) 

where m w max denotes the largest pion mass at which the expansion in used in a 
fit. A potential convergence criteria is to ensure that 
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a 4 m 7r,max 



(8) 



remains small. Here, Eq. ^ is assumed to be working to just leading nonanalytic 
order, ie. fully renomalised to 0(m^). The mt term is introduced in the fit to 
reduce the sensitivity to the ultra-violet cutoflP. It is desired that this term only 
contribute of the order 1% of the total result — corresponding to roughly 10% of 
the term, suggesting an initial test scale of TZ < 0.1. 



4. Rho-meson in quenched-QCD 

Quenched QCD offers a useful framework to test methods and techniques where the 
cost of dynamical simulations remain too computationally expensive. Further, in 
the case of the rho-meson, neglecting the dynamical sea quarks removes the decay 
channel of the rho. Thereby, one does not have to deal with complications arising 
through extrapolating through the 2-tt decay threshold^. 

On the assumption that disconnected contributions in the uj propagator are 
negligible, the only leading one-loop diagrams contributing to the p mass are those 
associated with the flavor-singlet rf . The rj' only appears in the low-energy EFT as 
a result of quenching, whereby it remains degenerate with the Goldstone pion. 

Our expansion for the rho-meson mass in quenched QCD is summarized as 

m p = ^m 2 p + S (9) 

where the residual, unrenormalised expansion is defined by 

rhp = a + a 2 m^ + a 4 m^ (10) 

and the chiral loop corr ecti ons are denoted by S. Further details will appear in a 
forth-coming manuscript!^. 



5. Preliminary investigation 

Here we wish to test the hypothesis to determine A by constraining the quantity 
TZ, Eq. ((8]). To do so, we construct some test, or pseudo, data in order to perform 
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hypothetical fits. Firstly, we fit Eq. ^ to new, unpublished, quenched, overlap 
results of the Kentucky GroupO^ using a dipole regulator of scale A = 0.8 GeV. 
This input simply enables a scale to be set for the test data. A more detailed study 
would investigate alternative constraints at this initial step. 

With this fit we construst a dense set of pseudo-data between the physical pion 

mass and SOITL6 UppGr value T2<7r,max- 

We then refit this pseudo-data with alternative 
regularization scales and assess the convergence criteria. The convergence test is 
plotted for varying A for a test curve which using an upper m ffimax = 0.3 GeV. 
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Fig. 1. Convergence criteria, 1Z, plotted against regularisation scale. 

The first indication of support for the procedure is the failure to produce a 
converging series if the regularisation scale is chosen as too small. This reflects the 
fact that the analytic expansion cannot describe the analyticitics of the data. The 
existence of a lower bound as determined by this criteria is then a signature that 
the chiral nonanalyticities are present in the data. 

There is a window of A values where the series is satisfactorily converging to meet 
the demands of 1% precision. This happens to be in the vicinity of the originally 
input scale of A = 0.8 GeV and success of the fit should come as no surprise. From 
our initial tests using pseudo-data constrained using A = 2.0 GeV, this convergence 
criteria reflects this and chooses an optimal regulator in this vicinity. Thereby, with 
good enough data, the intrinsic preferred scale can reveal itself through this criteria. 

The curve indicates that should a large regularisation scale be attempted to fit 
our pseudo-data, the residual series is not so well converging. Even here, with an 
upper pion mass of 300 MeV, truncation of the residual series suggests a (0.25) 2 or 
6% uncertainty. 

Another desirable feature of this procedure is that as m IjmH is reduced, the 
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sensitivity to A must disappear through the additional factor of max in 1Z. 
Thereby, within the PGR the lower bound will remain, yet no significant upper 
bound will persist. 

On face value, this looks like a reasonable procedure to allow lattice data itself to 
make an ab initio determination of an optimal FRR scale. A weakness is apparent in 
considering the renormalisation of the residual series. For this argument, considering 
just the renormalisation of the single-hairpin graph. The Taylor expansion of this 
diagram, up to normalisation, can be expressed as^ 

b + b 2 ml + ml + bimi + ... (11) 

The leading term bo behaves as A 3 , while £>2 ~ A and 64 ~ A -1 . Since for any fit the 
renormalised parameters are essentially stabilised by the data, this indicates that 
for very large A the residual series coefficients will behave as a{o,2,4} ~ 
Thereby, a-i will diverge and 04 saturate to a constant, and consequently our de- 
scribed convergence criteria 1Z will approach zero, regardless of the data. 

While a promising approach, it appears that the criteria described by Eq. (JSJ) 
does not lead to a conprehensive test of convergence. A potential modification 
could be to normalise the a^m% max term to a renormalised quantity, such as the 
rho mass. We anticipate further studies in this direction will lead to a reliable scale 
determination procedure that will facilitate ab initio studies with FRR for all A. 
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